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Abstract 

In energetic heavy ion collisions, if quark-gluon plasma is formed, 
its hadronization may lead to observable critical fluctuations, i.e., 
DCC formation. The strength and observability of these fluctu¬ 
ations depend on the initial state. Here we study the canonical 
ensemble of initial states of chiral fluctuations in heavy ion colli¬ 
sions and the probability to obtain observable domains of chiral 
condensates. 

PACS number(s): 12.39.Fe, 12.38.Mh, 25.75.-q 


I. INTRODUCTION 

In a recent article the initial states of Disoriented Chiral 
Condensate (DCC) were studied, the importance of non¬ 
vanishing chiral ’angular momentum’ was pointed out, 
and the microcanonical ensemble of initial states were 
evaluated for fixed energy and chiral angular momentum 
fli. These results indicated that the observability of DCC 
formation is a delicate problem, and it strongly depends 
on the reaction mechanism of heavy ion collisions. 

Since at this moment reaction models which would pre¬ 
dict a chiral angle distribution at the onset of hadroniza¬ 
tion do not exist, we assume that thermal equilibrium 
is reached in the supercooled plasma at the freeze-out, 
Tpo = Tq < Tc, which coincides with the beginning 
of the process of hadronization Consequently we 

also assume that at this moment the chiral degrees of 


freedom in the system are thermalized and thus can be 
characterized by a canonical distribution corresponding 
to the freeze-out temperature. To, and to the volume of 
the system, V. (This latter quantity is vital because we 
study fluctuations in a small, finite system!) Unless fur¬ 
ther constraints on the chiral degrees of freedom can be 
established from more detailed reaction models the gen¬ 
erated canonical distribution of chiral initial states will 
represent the statistical distribution in heavy ion exper¬ 
iments. 

Having obtained this distribution we will analyze which 
initial configurations lead to strongly fluctuating domains 
of DCC, and what are the overall chances to detect these 
experimentally. 

For our studies, just like ref. 0, we have chosen the lin¬ 
ear sigma model because it contains the relevant chiral 
degrees of freedom and it is still relatively simple. 

II. THE LINEAR SIGMA MODEL 

The linear sigma model involves the isodoublet of the 
light quark fields, 4*, and a real, four-component chi¬ 
ral field, $ = (cr, if), where tr is a scalar field, while 
TT = (tti, 712, TTa) IS a triplet of pseudoscalar pions. The 
Lagrangian density is [| 

£ = ^ — g (a + 775 ^ 77 )] 47 

+ I (di^cr) (d'^cr) + I (a^Tf) (a'^Tf) - U (cr, tt) , (1) 

with the so-called ” Mexican hat” potential energy density 

u {cr, tt) = ^ -b TT^ - -b - Ha + Uo{T) , 

( 2 ) 

where the last term, Uq{T), ensures that the absolute 
minimum of the potential energy density is 0. Without 
the Ha term, this Lagrangian density is invariant under 
the chiral SUl{‘2) <S) SUii{2) transformations. 

The parameters in this Lagrangian are chosen in such a 
way m that in normal vacuum at T = 0 chiral symmetry 
is spontaneously broken and the expectation values of the 
meson fields are 

(ct) = /,r, (tt) = 0 , (3) 

where /,r=93 MeV is the pion decay constant. In the 
vacuum with broken chiral symmetry pions represent soft 
“azimuthal” excitations of the chiral field. To have the 
correct pion mass in vacuum (T = 0), m,r=138 MeV, one 
should take 


*We shall use ?i = c = k = 1 units. 
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H = Uml . 


v^ = n- 


A 


( 4 ) 


The parameter A is related to the sigma mass, = 
2A/^ + TO^, which can be chosen to be about 600MeV 
(then A « 20). Sigmas represent stiff, “radial” excita¬ 
tions of the chiral field. The remaining coupling con¬ 
stant g can be fixed by the requirement that the effec¬ 
tive quark mass, + '^^), in broken vacuum 

coincides with the constituent quark mass in hadrons, 
which is about 1/3 of the nucleon mass, mjv- This gives 
g « ijUN/'i)/f-K = 3.36. 


tensor has the form § Tik = {e + p) UiUk -p gik, where 
e is the energy density, p is the pressure in the co-moving 
frame, and Ui is the flow velocity of the matter. Thus, in 
the LR frame the energy momentum tensor is diagonal, 
i.e., 


Tj; = diag (e, -p, -p, -p). 

The conservation laws lead to 

dep-\-e dp dp dp 

dT~^ T ’ dx dy drj 


( 8 ) 

(9) 


III. IDEAL HYDRODYNAMICS IN THE LOCAL 
REST FRAME 


These, together with the charge conservations 

A {PcU") = 0 , 


When solving the equations of motion corresponding to 
our Lagrangian we shall assume a one-dimensional 
Bjorken scenario^]. Consequently, we shall only con¬ 
sider field evolutions having one-dimensional scaling 
property and we shall work in the corresponding local 
rest (LR) frame. Therefore, let us briefly review here the 
properties of a scaling hydrodynamical expansion in the 
LR frame. 


A. Framework: the Bjorken model 


The local rest frame is a curvilinear coordinate system 
where the matter is at rest. For a longitudinal scaling 
expansion along the z-axis the LR frame coordinates, 
= (r, x,j/,? 7 ), are given by the c.m. coordinates, 
= {t,x,y,z), as 


T = 


and 



t + z 
t — z 


( 5 ) 


Thus, the invariant space-time interval in the LR frame 
is ds = dr^ — dx'^ — dy'^ — r^dg^, therefore the covariant 
metric tensor is 


g.fe = diag(l,-l,-l,-r^). (6) 


This metric tensor gives rise to three non-vanishing 
Christoffel symbols: Tgg = r, Tgg = rj^g = 1/r. 

The equations of ideal hydrodynamics are the conser¬ 
vation laws given by the requirement that the energy 
momentum tensor, T*, has a vanishing divergence, i.e., 
i Expressing this divergence in terms of the 
partial derivatives and the Christoffel symbols, we arrive 
to the following four equations: 


rjiK _ 


tL = 


dT^ 

dx^ 

dT^ 

dx^ 






f)rpk 1 

rpk _ 1 I rpO 

Tk _ ^"^3 rp3 /y\ 


In local thermal equilibrium, if we neglect the effects of 
viscosity and heat conductivity, the energy momentum 


(where c runs for all conserved charges, and Di stands 
for the covariant derivative) and the equation of state 
(EoS) which connects the pressure, energy density, and 
the proper charge densities pc [pc = dNc/{dx dy dg)\, 
form a closed set of equations, which can be solved for 
given initial and boundary conditions. In the LR frame 
the charge conservations simplify to 



( 10 ) 


The equation of the Bjorken model can be obtained by 
imposing the boost invariance and transverse transla¬ 
tional invariance of e, p and Pc- Then, by Eqs. (10), the 
charges are constants of motion: Pc{t) = Pc('n 3 ). Eqs. 
(^) simplify to an initial value problem for an ordinary 
differential equation in r: 


de(r) p{t) + e(T) ^ ^ 


( 11 ) 


dr T 

with e(To) = eg, where p = p{e, pco) is given by the EoS. 


B. Scaling expansion in the linear sigma model 

Let us consider the sigma model without quarks (5 = 0). 
Then the Lagrangian density (^ , in a general curvilinear 
coordinate system, is 

£ = ^(d,<l>r)g^Udk<^r)-U(<Pr,T), ( 12 ) 

where di = ^ 7 , = (o',!?), there is an implied summa¬ 

tion for indices occurring twice (i, k and r), and 

U{a, 7?, T) = ^ + 7 f^-v^ + - Ha + Ug, 

where a:® is a space-time coordinate a;® = (r, a;, y, ry). 
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The energy-momentum tensor, Tik, can be obtained from 
the Lagrangian density using the relation^] H] 

_1 /_ d 

2 '^ ^ dxi dgi’^ ’ 

dx^ 

where g stands for the determinant of the metric tensor, 
gik- A straightforward calculation gives 

T^k = {di^r){dk^r)-Cg^k- ( 13 ) 

This is, unlike the energy-momentum tensor in the 
Bjorken model, generally not diagonal. However, if we 
impose boost invariance, the off-diagonal elements in the 
energy-momentum tensor vanish. From Eqs. (|l^ ) 

and (g), we can easily obtain the diagonal elements: 

Tq = — (c^T'hr) + U{^r) = Skin + Spot = S 

— 2 (^rdlr) (^*7.) = Skin Spot = P 

Ti = Ti=Tl, ( 14 ) 

where Ckin = 5 (9r$r) (clr^’r) is the kinetic energy den¬ 
sity and Cpot = U is the potential energy density, and we 
introduced the total energy density, e = Skin + Spot, and 
p = Cpot — Skin- These are functions of the proper time r 
as the fields follow the equations of motion (EoM). 
Obviously, the energy momentum tensor now has the 
form (^, therefore the conservation laws (^) and all their 
consequences are also valid now. Particularly, this leads 
to 


de _ S p _ 2 Skin 

dr T T 


(15) 


which can be seen as a generalization of the usual en¬ 
ergy conservation for our expanding system. We used 
this equation as a check of the numerical solution of the 
equations of motion. 

To obtain how the total energy evolves, we have to take 
the expansion into account. The surface element of the 
T = const, hyperplane is dT, = rdpdxdy. Therefore, the 
volume of our longitudinally expanding system is propor¬ 
tional to T, i.e., dV/dr = VjT. Now, with the help of the 
conservation laws ([III), the time derivative of the total 
energy, E = Ve, can be expressed as 


dE Epo^ ddkin dV 

dr T ^ dr ^ 


(16) 


where Ekin = V Ckin is the kinetic energy and Epot = 
V Cpot is the potential energy of our expanding system. 


^Our sign convention in the space-time interval is different 
from that of ref. 0; therefore, we have a minus sign on the 
left hand side. 


The interpretation of Eq. (^^ could be similar to that 
of the corresponding equation for the Bjorken model, Eq. 
&■ However, for the sigma model we can easily have 
a negative pressure, which corresponds to an increase of 
the total energy duriM expansion! This is the case if 
Epot > Ekin in Eq- (|l6|), which is typical for the early 
evolution in the quench scenario. We also note that, as 
pointed out already in ref. [Q, Eq. ( |l^ means that we 
have an isentropic (adiabatic) expansion. 


IV. EQUATIONS OF MOTION 


In the mean-field approximation, ignoring all-loop con¬ 
tributions and considering a and if as classical fields, we 
can obtain the following equations motion : 

dpd^(j{x) -I- A [cr^(a;)-|-7f^(a;)—r:'^(a;)] cr(a:) 

-H = -gps{x), 

dfj,d^Tr{x) + A [cr^(a;)-|-7f^(a;)--y'^(a;)] 7f(a:) = -gpp{x), 

(17) 

where v'^{x) = — T'^{x)j2. Here ps = {qq) and 

Pp = i{qj 5 Tq) are the scalar and pseudoscalar quark den¬ 
sities, which should be determined self-consistently from 
the motion of q and q in the background meson fields. 
The scalar and pseudoscalar densities can be represented 
as @ 


Ps{x) = ga{x)a{x), pp{x) = ga{x)Tf{x), (18) 


where a{x) is expressed in terms of the momentum dis¬ 
tribution of quarks and antiquarks f{x,p): 


a{x) = Vq 


d'^p 

{2'xTiY 


‘2S{p^Pp 


'sn?{x))f{x,p) 



d?p 

E{x,p>) 


[nq{x,p>) -f nq(a::,p)] . 


(19) 


Here Vq is the degeneracy factor of quarks which equals 
to 12 in our case (flavour x spin x colour). 

Let us consider a one-dimensional Bjorken scaling and 
assume that the quark distribution is an ideal Fermi dis¬ 
tribution at the freeze-out time tq (following [||J|] we can 
take To ~ 7fm/c) while the subsequent distribution for 
time T > Tq is given by the solution of the collisionless 
Vlasov equation |^. Then the source term, a{x), can be 
cast into 


a(x) = a{T) = 


X dq q 


VnOi 


TT^Vl ~ oi^ 


rjiZ 

-L n 


arcsin^(1 - a^) / +” 

exp (^^/q^ + m^{To)/T^ - 


+ 1 


where a = tq/t, Tq and ptQ are the freeze-out tempera¬ 
ture and chemical potential. For the po = 0, mq = 0 case 
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this expression can be integrated analytically , while in 
the general niq 7 ^ 0 , /tq 7 ^ 0 case it can be evaluated nu¬ 
merically using an eight-point Gauss-Laguerre integral 
formula [^. This way we can follow the solution of the 
EoM numerically from an arbitrary initial condition at tq 
during the whole symmetry breaking process. For sim¬ 
plicity, we shall assume that quarks and antiquarks are 
in chemical equilibrium, i.e., we take /xg = 0 . 

To solve the EoM, we need to know the temperature field, 
T{x). However, following ||l^,|^ we may assume that the 
cooling is very fast and then we may take T(x) = 0. 
This is the so-called quench scenario. Certainly, the 
quench is an idealization since in a real collision ther¬ 
mal fluctuations will not cease instantaneously. Now 
the EoM has no reference to any temperature but the 
initial condition, which we choose to be at the freeze- 
out of the quarks, must represent the freeze-out tem¬ 
perature, To « 130 — 140MeV. We emphasize that our 
scenario differs from the original quench because, even 
though the thermal fluctuations of the fields are ne¬ 
glected, we still have an interaction with the expanding 
quark/antiquark background. One could choose other 
possible non-equilibrium scenarios, for example, anneal¬ 
ing 

Before the quench, the system can typically be found 
somewhere around the minimum of the potential energy 
density. However, when the quench occurs, the position 
of this minimum suddenly changes to the physical vac¬ 
uum, (cr, 7 ?) = (/^, 0). Thus, the system finds itself on the 
slope of the ’’valley” of the potential, rolls down towards 
its new minimum, and oscillates around it . 

The oscillations are damped because, due to the expan¬ 
sion, there is a first order derivative in the EoM (0). 
therefore the fields finally relax to the physical vacuum, 

(A.o)- 

Studies show that the main effect of the quarks 

is an additional damping; the oscillations of the chiral 
fields are reduced by about 40 — 50%. This is due to the 
fact that during the transition a portion of the energy 
of the system is used up for the increasing quark mass. 
Also, the phase transition happens about 0.5fm/c slower; 
it takes about 2.5fm/c. The presence of quarks has a 
strong influence on the onset of instabilities and domain 
formation in the initially homogeneous system. This we 
shall address in a separate paper. 


V. THE QUESTION OF THE INITIAL 
CONDITION 


The physical choice of initial condition or ensembles of 
initial states was not studied in detail in earlier works, 
although this is a vital problem if we want to predict the 
observability of DCC fluctuations. An early study 0 
followed by refs. | W |] assumed the initial condition: 


CTo =0, TTo = 0, 

do = IMeVc/fm, ttq = (5, 0,0)MeV c/fm, 


with To = Ifm/c. This choice of the chiral field configu¬ 
ration corresponds to the chirally symmetric phase rep¬ 
resenting very hot quark-gluon plasma, i.e., when H = 0 
and T > Tc = \f2v (~ 124MeV). The initial value of 
the fields was assumed to be the value of the condensate, 
which is zero in this case, i.e., (cr) = 0 and (7?) = 0. 
However, to get the correct value of the pion mass in 
the physical vacuum, we do need the explicit symmetry 
breaking term Ha. Hence, in the minimum of the poten¬ 
tial the condensate is never zero, though it gets close to 
zero for temperatures above 250 — 300MeV (see Fig. |^). 
However, at the freeze-out temperature, Tq = 130MeV, 
we have a large condensate of 40MeV. 

Another important factor is that our system is finite, 
therefore thermal fluctuations are of real importance. 
Due to the fluctuations, the system is at the minimum 
of the potential energy density only on the average, and 
the probability for finding it somewhere around this min¬ 
imum can be considerable, especially for very small sys¬ 
tems and/or high temperatures. Therefore, the time evo¬ 
lution starting from the minimum of the potential energy 
density is not fully representative for the majority of the 
realistic initial states. 

Furthermore, the system has a considerable kinetic en¬ 
ergy at the temperatures under consideration. This 
should be reflected by the choice of the initial value of 
the field derivatives. Our analysis will show that the ini¬ 
tial condition has too little kinetic energy for the 
temperatures we are considering, even if we go down to 
as low as T = 130MeV. Thus, the picture of a system 
on the ’’top” of the potential started with a little ’’kick” 
cannot be realistic. 

These considerations lead us to propose an initial condi¬ 
tion, which is closer to the physical picture we have at 
the time of the freeze-out. Before the freeze-out the sys¬ 
tem can be approximated to be in thermal equilibrium 
at a given temperature, T, and it has a given volume, V. 
If so, then the initial condition should be chosen accord¬ 
ing to a canonical thermal ensemble, i.e., the probabil¬ 
ity p of a configuration {a, 7?, a, 7?) is proportional to the 
Boltzmann-factor, P: 


p{a, TT, a, tt) (X exp 

= P(ct, TT, CT, 7?) 


V ■ ' 

-—H{a,Tt,d,Tt) 


( 21 ) 


where TL is the Hamiltonian density, the dot stands for 
djdr, and we assumed that the system is homogeneous 
and we have taken the Bjorken scaling granted. The 
Hamiltonian density consists of three terms, TLkin, 'Hpot 
and Tiq, which are respectively the kinetic and potential 
energy density of the chiral fields and the energy density 
of the quarks (with p = 0): 

'Hq = Vq J + fq{x,]^] 


( 20 ) 
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‘iVq 

{2^f 


(fp 


exp 


\/p^ + m? 
+ rn? jT 


+ 1 


ni • 2 

^kin — 2 ^ 2 ^ 

'Hpot = U{a,n) = -^ (cr^ - v''^)^ - Ha + Uo , 


( 22 ) 


where the effective quark (antiquark) mass is given as 
TO^ = g^{a^ + TT^). 

Similar approaches can be found in the literature where 
Gaussian distributions were considered [llT|-H- How¬ 
ever, as the potential energy density is not quadratic and 
we also have the quark term, the distribution of the chiral 
fields is clearly non-Gaussian. We shall give a comparison 
to the initial condition of refs. 

{ 4 >) = (^) = 0, (02) = ^;2/4, (02) = ^;2/lfni2. (23) 

A forerunner of the present analysis is ref. 0 where we 
considered initial conditions with large chiral ’’angular 
momenta”, /, i.e., states where the chiral fields circle 
around the minimum of the potential energy density hav¬ 
ing a trajectory which is far from the minimum. As we 
have shown there, these states as well lead to DGG for¬ 
mation during the fast linear expansion. However, the 
magnitude of I remained an open question. Now we shall 
extend the analysis by including the initial distribution 
of finite chiral ’angular momentum’ states. 

There are several methods to generate the ensemble of 
the initial conditions. In ref. ^ we needed a microcanon- 
ical ensemble around a specified value of the ’’angular 
momentum”, I. This ensemble was obtained by sam¬ 
pling the time evolution of the fields evenly in the case 
with explicit symmetry breaking but without quarks and 
expansion (Fig. ||). In this case I is not a constant of 
motion but oscillates in a certain interval. This method 
is based on the assumption of the ergodicity of the sys¬ 
tem, i.e., that during the time evolution the system gets 
arbitrarily close to any point in the available phase space. 
It is also possible to generate a canonical ensemble based 
on ergodicity. In ref. a Langevin type equation of 
motion was used to generate such initial states, where a 
connection to an external heatbath of a given temper¬ 
ature was assumed. Then, based on the ergodicity, the 
states at different times correspond to a canonical ensem¬ 
ble of the system. 

In this work we generate our initial canonical ensem¬ 
ble using a different method, namely the Monte-Carlo 
method. Consequently, we shall use a different formal¬ 
ism from that in ref. or ref. ||^. 


A. The generation of the canonical distribution 

To generate the canonical distribution d^), we used a 
simple generalization of the so-called rejection method 
0 for a multidimensional probability distribution on a 


boundless manifold. A convenience of our method is that 
it does not require proper normalization of the probabil¬ 
ity distribution, i.e., the evaluation of the normalization 
factor. 

To apply the method we need to find a function, 
F(tT, if, (T, if), which is everywhere above P{a, 7f, ct, tI) and 
is proportional to a distribution, /(cr, if, ct, 7f), which we 
are able to generate. It is simple to prove that if one 
generates a (8-dimensional) random deviate, (cr, if, tr, if), 
of / and a uniform deviate y between 0 and F’(cr, if, d, if), 
then the (tr, if, d, if, y) points will be uniform in the (9- 
dimensional) volume under the graph of F. If now one 
keeps the point (u, if, d, if, y) whenever it is under P and 
rejects it if not, the accepted (u, if, d, if) points will have 
the desired p{a, if, d, if) distribution. The fraction of ac¬ 
cepted points, i.e., the efficiency of the method is 

f • ■ • f P(o', if, d, if) da da d^ir 

^ -:- 

I ■ ■ ■ S 

To make it easy to generate the / distribution corre¬ 
sponding to F, we shall choose F to be proportional to 
a product of Gaussian deviates: 

F{a,iT,a,iT) = 

= Fi(cr,if) X F 2 (d,if). (24) 

Now / will simply factorize into one-dimensional Gaus¬ 
sian distributions of the chiral fields and their derivatives. 
These Gaussian deviates can be generated using the rou¬ 
tine gasdev described in ref. 0- The random number 
generator we used is the routine ran2 in ref. ||^ . 


1. The choice of the eomparison function 

We have to choose the parameters, C, Ci, C 2 and in 
F such that 


F{a, if, d, if) < P{a, ff, d, if) (25) 


for all values of the chiral fields and their derivatives. 
The coefficient C 2 is determined by the kinetic energy 
density term because the field derivatives appear only 
there. Since it is a quadratic function of the derivatives, 
it is straightforward to choose C 2 = VjFr, i.e.. 


F 2 (d, 7f) = exp 



(26) 


The fields are in the other two terms of the Hamiltonian 
density, which are not quadratic, therefore we need ap¬ 
proximations. With the help of Eq. (p^, inequality ( [^ ) 
can be written as 

V 9 0 

y (Idpot + T-iq) >Ci{a- am) + Giff -|- const., 
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thus we have to approximate the potential and quark 
energy density terms with a quadratic expression from 
below. 

The quark energy density contains the fields in a rather 
intractable integral. It is a monotone decreasing func¬ 
tion of the quark mass and it is in the range between 0 
and 257MeV/fm^. It turns out to be sufficient to take a 
simple, crude lower bound of the form L — Mtv?, where 
M >0. 

It is enough to deal with the quartic term in the poten¬ 
tial energy density (p^, since the remaining Ha term 
can easily be combined with any quadratic expression of 
a into a final Ci {a — amf' quadratic expression. The 
quartic term has the general form C (a^ ± l) , where 
0 ? = (cr^ -I- TT^) /|u'^|, the upper sign (-I-) corresponds to 
T > Tc while the lower one (—) to T < Tc, and C does 
not contain the fields. We approximate this expression 
by a quadratic one using that (a^ ± l) > — k, if 

0 <l <2y/k + l±2. (27) 


Putting everything together, the first part of the com¬ 
parison function is 


Fi (a, if) = exp 


X exp 


j V X^l'\v'^\ 
1 ~T 4 

'V / H^ 

f 



2H y 


x'^F'^k 


(28) 


with /' = Z-4MP/A2 |u'2|. 

It is crucial to choose the parameters I, k, L and M such 
that we have a reasonable efficiency. Clearly, to get the 
highest efficiency we have to minimize the area under 
Fi. It is easy to see that for a fixed I we need k to be 
minimal, thus we have to take the upper limit in Ineq. 
(p7|), I = 2\/k F 1 ± 2. Since the quark energy density 
is a convex function of m^, for a given M the optimal, 
i.e., the largest L is given by the tangent line that has 
a slope (—M). Now we have only k and M left. For a 
given M, the minimum condition for the area under Fi 
leads to a third-order algebraic equation for k that we 
solved numerically. M we set by trial and error. Having 
the above optimization we obtained an efficiency of 30— 
60%. 

The quark energy can be expressed as 


Hq 


CO 




fe-i 





where stands for the modified Bessel function of the 
second kind of n-th order. Certainly, we cannot take 
into account all terms here but that is unnecessary as 
the sum is convergent. We summed up the terms until 
the relative change of the sum became less than I 0 “® 
(this corresponds to 4 — 16 terms). 

The above formula cannot be used directly for m = 0 and 
for T = 0. For m = 0 we can obtain the expression |^] 


Hq = I'qT'^/120 . For T = 0 the inte^and and thus 
the quark energy density is zero (see ref. |^, Section 7.3). 


B. An advantage of the symmetries of the 
Hamiltonian 


Though it does not take long to generate the ensemble of 
the proper initial conditions, it is very time-consuming 
to compute the time evolution from the several thousand 
initial states. However, it is possible to exploit the sym¬ 
metries of the Hamiltonian to reduce the computation. 

It is easy to verify the following symmetry: if a trajectory 
(ct (t) , tt (t)) satisfies the EoM (p^, so does the trajec¬ 
tory (cr (r) ,07? (t)), where O is an arbitrary, 3x3 ro¬ 
tation matrix in isospin space. Obviously, the rotation 
implies that the initial condition as well should be ro¬ 
tated. Now, if we have data for a given initial condition, 
we at once know the results for all other initial conditions 
obtained from this given one by rotations; we simply have 
to rotate the pion fields using the same rotations. 

In this sense all rotationally equivalent initial conditions 
are the same. We may use this to reduce the number of 
parameters in the initial condition. Originally these are 
8 (four fields -I- four derivatives). However, taking into 
account the rotations we may fix two pion fields and the 
sign of the third one, e.g., we can take tti = 7 r2 = 0 and 
TTa > 0 , and we can also fix one derivative, e.g., 1(2 = 0 . 
Thus, we have only 5 parameters left. 

Since our canonical distribution contains the pion fields 
only through and tt , it is convenient to reexpress the 
above parametrisation using the length of the isovector, 
| 7 r|, the length of the derivative, | 7 f|, and the angle, 9, 
between 7 ? and tt. Their relation to the previously used 
parameters is trivial: 773 = |7r|, ( 7 r 3 , 7 ri) = |7r|(cos0,sin0). 
Now we may realize the distribution in spherical coordi¬ 
nates: 


exp 

exp 


V 

-—n(a,Tf,a,a) 

|5?|,d-, |ir|) 


da dif da dif = 


if’if da d\if \ da d\if\ 


X sin 0 dQ d^ sin d dd dip , 


(29) 


where the angles d and p correspond^ to if, while the 
other two angles, 0 and 4), correspond to 7?. The ori¬ 
entation of the two spherical frames can be chosen arbi¬ 
trarily. We take the choice when if is along the z-axis in 
the frame of if because then 0 is the angle between these 
two vectors. Exploiting the rotational invariance, the an¬ 
gles p, d and 4) can be fixed, e.g., to tp = -d = 4> = 0 , 


*Here, we mean the usual parametrisation tti = 
| 7 f| sin 19 cosp, 772 = | 7 f| sin 19 sin(p, and 773 = | 77 | cosi9, with 
p e [ 0 , 277) and i9 e [ 0 ,77); and a similar one for the vector tt. 
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and the 5 parameters needed for the initial condition are 
tj, |7r|, (T, |7f| and 0. We will only need the distribu¬ 
tion of the three angles, (p, d and d*, to reconstruct the 
complete, 8-variable distributions of any physical quan¬ 
tity from the results for the ensemble of the 5-variable 
initial conditions. 

The distribution d^ ) factorizes into three terms: a term 
depending on the fields and field derivatives, cr, d, |7r| and 
|7r|, a term depending on 0, and a third term containing 
the other three angles, $, d and (p. We may generate the 
first term using the rejection method introduced in the 
previous subsection. The second term can be obtained 
by transformation from a uniform deviate (see ref. [p^ ). 
This way, we can have a representative ensemble of the 
5 parameters needed for the initial conditions. We can 
easily obtain the ensemble for the remaining three angles, 
for example, using the rejection method again. By form¬ 
ing pairs from the elements of the two ensembles, (let us 
define this operation as the product of the two ensem¬ 
bles) we have a representative ensemble of the original 
canonical distribution. However, we need to run the cal¬ 
culations only for the first ensemble in the product, while 
the second one we apply on the results as rotations. This 
enables us to choose the second ensemble as large as we 
wish. Furthermore, there are parameters, for example 
the neutral pion ratio, that we can calculate analytically 
by taking into account the distribution function of the 
second ensemble! 

Let us show this in more detail. The joint distribution 
of the angles p and "d is independent of the other vari¬ 
ables in (|^), moreover, it is the uniform distribution on 
the three-sphere, given by siniD (I'd dp. To reconstruct 
the distribution of the tt field, we have to perform the 
rotations given by D and p on the states in the ensemble, 
which leads to a spherically symmetric tt distribution at 
all proper times. Therefore, even though unlike ref. 0 
we do not have an 0(4) isospin symmetry, for the neutral 
pion ratio, R, we obtain exactly the same distribution: 


V{R) = 


1 

2/r' 


C. Results for the ensemble of initial conditions 

In this subsection we shall present the results of the 
generated ensembles of initial conditions. These ensem¬ 
bles depend on the freeze-out temperature and chem¬ 
ical potential. To and po, and on the volume of the 
DCC domain. The former two are determined since we 
start from a 20% supercooled baryon-free plasma, i.e.. 
To = ISOMeV and po = 0. 

It is known |j^ that dX T k. T^. the correlation length 
in the linear sigma model is comparable to the pion size. 
To simulate these most probable initial domain sizes, we 
have chosen the value Vq = 0.5fm^. 

To give a comparison to the Gaussian initial condition 
(p^), we also tried Vq = 3.3fm^ since that gives exactly 


the same initial field derivative distributions. 
Unfortunately, our model does not include the dynamical 
evolution of the domain sizes. It is not clear what the 
survival rates for domains of different sizes are and it 
may occur that relatively large domains have the best 
chance to survive. That is why we chose a third value, 
Vo = lOfm^, for our analysis. 

All three initial ensembles contained 15,000 states. 


1. Various distributions of the initial ensemble 

a. The distribution of the chiral fields. Fig. ^ shows the 
joint distribution function]^ of the chiral field a and the 
length of the isovector tt in the initial ensemble. As we 
expected, the peak of the distribution is not at cr = 0, 
however it is not at the minimum (cr « 40MeV) of the 
potential energy density, either. This is because of the 
quarks. The quark energy density is a monotone decreas¬ 
ing function of cr^ -I- rf^, therefore it is more favourable 
energetically to have a larger condensate, « 86MeV. 

For bigger volumes the distributions are sharper. The 
reason is that fluctuations are smaller in larger systems. 
Thus, the probability for the system to be far from the 
minimum of the potential energy density is larger for 
small systems. This increases the probability of having 
large oscillations, which helps the observability of the do¬ 
main. However, this increase may be compensated by the 
fact that small domains emit fewer pions. Unfortunately, 
our one-dimensional model does not allow us to decide 
this question because it includes only a scaling expansion 
of the domains but cannot account for domain evolution 
due to instabilities. 

As a comparison. Fig. ^ shows the cr— |7f| distribution for 
the microcanonical ensemble shown in Fig. 7 in ref. [^. 
The states in the ensemble are concentrated in a circular 
band, which is rather different from our peaked distri¬ 
butions. This is a consequence of the high initial chiral 
’’angular momentum” and the constraint that these ini¬ 
tial states all have the same energy. 

b. The distribution of the quark/antiquark effective mass. 
The initial distribution of the chiral fields implies that 
the initial quark mass is very unlikely to be zero, which 


^ These distributions were obtained by introducing a rectan¬ 
gular grid in the a — Ivf | plane and counting the states in each 
grid cell with a weight. The value of the distribution 
function in a given grid cell was approximated by the ratio of 
the total for the cell to the area of the cell. Then we divided 
by the total number of states to normalize the distributions to 
unity. Since the I/tt^ Jacobian enhances statistical fluctua¬ 
tions for small |7f| values tremendously, we had to take a huge 
ensemble of 600,000 states for these plots. These fluctuations 
also depend on the cellsize; the distributions are smoother for 
bigger cells. We have chosen a cellsize of 20MeV in both the 
a and I if I direction. 
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again questions the appropriateness of the choice (|20[). 
The actual distribution^ of the effective quark mass is 
shown in Fig. As the initial volume decreases, the 
distribution broadens and the average quark/antiquark 
mass increases. 

c. The distribution of the proper time derivatives of the 
fields. By Eqs. ( ^l|) and (|^), the proper time deriva¬ 
tives of the chiral fields follow Gaussian distributions 
with the same width, y^To/Vo- For a given freeze-out 
temperature, the width is inversely proportional to ^/Vq. 
In ref. [|^ and followers | P^ , P^JT8 Hi^ this width is as¬ 
sumed to be order of v, interpreted as u/lfm, which 
is 88MeV/fm. This is the same order of magnitude as 
our results and would correspond to Vq = 3.3fm^ at 
To = 130MeV. To obtain a significantly lower width, 
one has to take an extremely low freeze-out temperature 
and/or a very large domain size. 

d. The distribution of the quark, potential and kinetic en¬ 
ergy densities. The distribution of the quark, potential 
and kinetic energy densities of the initial conditions can 
be seen in Figs. |^, and |^. We can see that the kinetic, 
potential and quark energy densities have the same order 
of magnitude. 

Due to the increasing fluctuations, all distributions 
spread as the initial volume decreases. Furthermore, the 
peaks of the kinetic and potential energy density distribu¬ 
tions shift to higher values. This is in agreement with our 
picture: more fluctuations result in more kinetic energy 
density, and thus the system climbs higher and higher on 
the wall of the valley of the potential energy density. On 
the other hand, for the quarks the peak shifts to lower 
values for increasing initial volumes. This is due to the 
fact that fluctuations increase the quark/antiquark effec¬ 
tive mass. 

The distribution of the kinetic energy density ( ^| ) can be 

calculated analytically. Using that TCkin = ^ /2, from 
Eq. (p9|) it is not difficult to obtain that 


f {'Hkin) dTlkin = 


exp 


rr 

-^0 


Tik in dUk in 


A simple calculation yields that the maximum of the ini¬ 
tial kinetic energy distribution is at Ttkin = Tq/Vq. This 
value is 2 — 3 orders of magnitude larger than the kinetic 
energy density corresponding to the initial condition (20), 
which is about 0.066MeV/fm^. 


2. The evolution of the ensemble 


a. The evolution of the ensemble of the chiral fields. 
Figs. H and || show the time evolution of the ensem¬ 
ble average and the 70% fraction of the chiral fields, a 
and TTi, for freeze-out times tq = Ifm/c and 7fm/c. The 
other two pion fields, tt 2 and have a similar evolution 
to that of TTi because of the rotational symmetry. 

Here, by average we mean the arithmetic mean, while the 
”70% fraction” is the interval that is centered on the av¬ 
erage and contains 70% of the ensemble. One could also 
use the usual estimate of the standard deviation but, be¬ 
cause the distributions are non-Gaussian, that would not 
correspond to the same fraction of the ensemble at dif¬ 
ferent proper times (we determined that in our case one 
standard deviation means an at least 60% fraction, while 
two standard deviations correspond to at least 90%). 
The 70% fraction helps us pick out very unlucky events, 
i.e., those which are very far out of the interval. These 
intervals are narrower for bigger systems, due to the well- 
known fact that fluctuations are smaller in bigger systems 

. _ , 

The time evolution from the initial condition (P0| ) with 
To = Ifm/c is rather uncharacteristic for the majority, at 
least 90%, of the ensemble, as from time to time the a 
field is away from the average by about two standard de¬ 
viations. For To = 7fm/c the situation is slightly better - 
the above fraction is 80%. This is because the damping 
is reduced due to the slower expansion, therefore the dis¬ 
tribution becomes wider. We conclude, that the initial 
condition is representative only if the initial ensem¬ 
ble has larger fluctuations, i.e., we have a much higher 
initial temperature or smaller initial volume, which is 
physically unreasonable. 

However, the evolution of the pion fields is well in the 
70% fraction. Thus, the initial condition (po| ) does not 
overestimate pion oscillations. Furthermore, due to the 
small starting kinetic energy, the oscillations are deep in 
the 70% interval, showing that oscillations with larger 
amplitudes would be similarly probable. The oscillations 
are out of the 70% fraction only if the initial ensemble 
contains little fluctuations, i.e., if the initial temperature 
is very low and/or the initial volume is very large. 

Figs. ( p[) a nd also tell that the Gaussian initial con¬ 
dition ( p3|) results in a very similar ensemble evolution 
to that for Vq = 3.3fm^ and Tq = ISOMeV. The 70% 
intervals are almost identical, showing that both ensem¬ 
bles are concentrated in the same regions. The mean 
values are identically zero for the tti distribution as a 
consequence of the rotational symmetry. However, the a 
means are quite different because the distributions have 
different shapes and the initial distribution is peaked 
at cr = 0, which is not the correct value of the condensate 
at To = 130MeV. 


** All one-dimensional distribntions were calculated by di¬ 
viding the plotted interval into 40 — 50 equal sections and 
counting the number of states in each section. All these dis¬ 
tributions are normalized to unity. 







VI. CONCLUSIONS 


In this work we investigated chiral condensate evolution 
starting from a canonical ensemble of initial states. In the 
scenario we presented, the initial ensemble and its subse¬ 
quent evolution is characterized by three simple param¬ 
eters: the freeze-out temperature, Tq, the initial volume, 
Vqi and the freeze-out time tq. 

We compared our initial condition to the ’’roll-down” 
scenario, i.e., to starting from the unstable maximum 
of the potential with a little kick We found 

that the sigma field evolution in the ’’roll-down” scenario 
is unrealistic because the sigma field oscillates counter¬ 
phase, however pion oscillations are not overestimated if 
we choose reasonable values for the parameters, i.e., if 
To = ISOMeV, Uo - 0.5 - lOfm^, and tq - 1 - 7fm/c. 
Thus we confirm the suggestions and the findings 
[ p^ that observable DCC is a rare event in a canonical 
ensemble. 

We also made a comparison to the often-quoted Gaussian 
field and field derivative distributions with widths v and 
ri/lfm respectively 0 We found that these 

distributions are quite close to the choice Tq = 130MeV 
and Vo = 3.3fm^. However, this choice also accounts for 
the necessary non-zero mean of the a distribution corre¬ 
sponding to the non-zero condensate at a finite temper¬ 
ature, To ~ 130MeV. 

We believe that our simple scenario gives more insight 
on what the proper initial conditions could be. Though 
our results are similar to those obtained from other initial 
conditions, the ones mentioned above and their modifica¬ 
tions , we did contribute to reducing the arbitrariness 
in those. Our three parameters have clear physical mean¬ 
ing, and they give a simple way and a good handle on 
how to choose the initial field distributions. 

One should keep in mind that all our calculations were 
done in a one-dimensional Bjorken scaling picture. To ob¬ 
tain quantitative answers on domain sizes one inevitably 
has to take into account the spatial dimensions as well. 
Our test calculations based on the linear response method 
indicate that the fastest growth rates are about Ic/fm, 
thus the biggest domains are 3 — 4fm large. This is in 
agreement with earlier findings [|l^[l4[[D^ P,[l5| . However, 
we would like to emphasize that the only way to make real 
quantitative statements about the observability of DCCs 
is to go one step further and calculate observables, e.g., 
to obtain the spectra of pions emitted by the domains. 
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FIG. 1. The minimum of the potential energy density as a 
function of the temperature with our standard model parameters. 
Clearly, the condensate is non-zero even at high temperatures. 
The ’diamond’ indicates its value at the freeze-out temperature 
To = 130MeV. 



FIG. 2. The initial distribution of the chiral field cr and the 
magnitude of the isovector tt, for the microcanonical distribution 
in Fig. 7 in ref. [P. Unlike the peaked canonical distributions 
in Fig. the states are concentrated in a circular band. The 
initial ensemble contained 600,000 states and, to aid the eye, the 
distribution was complemented by a refiection on the |7r| = 0 plane. 





FIG. 3. The initial distribution of the chiral field cr and the mag¬ 
nitude of TT, for initial volumes Vq = 0.5fm^ (top), 3.3fm^ (middle) 
and lOfm^ (bottom), at the initial temperature Tq = 130MeV for 
an ensemble of 600,000 states. Manifestly, the distribution has its 
maximum at a non-zero g value and becomes narrower for a larger 
initial volume. To aid the eye, the distributions were complemented 
by a refiection on the |7f| = 0 plane. 
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FIG. 4. The initial distribution of the effective quark/antiquark 
mass for initial volumes Vq = 0.5fm^ (pluses), 3.3fm^ (crosses) 
and lOfm^ (boxes), at the initial temperature Tq = 130MeV for 
an ensemble of 15,000 states. The ’diamonds’ correspond to the 
Gaussian initial condition (^). As the initial volume decreases, 
mass fluctuations intensify and the average mass increases. 



FIG. 6. The initial distribution of the potential energy density 
( [ 2 ^ of the meson fields for initial volumes Vb = O.Sfm^ (pluses), 
3.3fm3 (crosses) and lOfm^ (boxes), at the initial temperature 
To = 130MeV for an ensemble of 15,000 states. The ’diamonds’ 
correspond to the Gaussian initial condition (^3|). As the initial 
volume decreases, the system is more and more probable to be on 
the walls of the valley of the potential energy density, far from the 
minimum. 



FIG. 5. The initial distribution of the quark energy density den¬ 
sity ( [ 2 ^ for initial volumes Vq = 0.5fm^ (pluses), 3.3fm3 (crosses) 
and lOfm^ (boxes), at the initial temperature To = 130MeV for 
an ensemble of 15,000 states. The ’diamonds’ correspond to the 
Gaussian initial condition ([ 2 ^. Since for smaller initial volumes 
the quark mass tends to be bigger, the quark energy density is 
then smaller. 
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FIG. 7. The initial distribution of the kinetic energy density 
of the meson fields for initial volumes Vb = 0.5fm^ (pluses), 
S.Sfrn^ (crosses) and lOfm^ (boxes), at the initial temperature 
To = 130MeV for an ensemble of 15,000 states. The distribu¬ 
tion for the Gaussian initial condition ( p^ ) is identical to taking 
Vb = 3.3fm^. The agreement with the theoretical curves is excel¬ 
lent. For smaller initial volumes fluctuations are larger, and thus 
the kinetic energy distribution is broader and the average kinetic 
energy density is larger. 
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FIG. 8. The time evolution of the ensemble average and ”70% 
interval” of the cr and tti fields for initial volumes Vb = 0.5fm^ 
(top), 3.3fm^ (middle) and lOfm^ (bottom), with a freeze-out time 
TO = Ifm/c, initial temperature To = 130MeV, for an ensemble of 
15,000. The shaded regions represent the 70% fraction centered on 
the mean value. The thin solid lines correspond to the evolution 
starting from the initial condition ([2o|). The evolution of the a field 
is rather uncharacteristic for the majority of the ensemble, showing 
that this initial condition is improbable. However, the oscillations 
of the pion fields are well within range, which is due to the rather 
small kinetic energy corresponding to (|20|). The middle figures also 
show the ensemble evolution for the Gaussian initial condition 
mean values — crosses, 70% interval - dashed. The evolution is close 
to that from Vb = 3.3fm^. 
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FIG. 9. The time evolution of the ensemble average and ”70% 
interval” of the cr and tti fields for initial volumes Vq = 0.5fm^ 
(top), 3.3fm^ (middle) and lOfm^ (bottom), with a freeze-out time 
ro = Ifm/c, initial temperature Tq = 130MeV, for an ensemble of 
15,000. The shaded regions represent the 70% fraction centered on 
the mean value. The thin solid lines correspond to the evolution 
starting from the initial condition ( |20[). The evolution of the a field 
starting from the initial condition (EOt) is still uncharacteristic for 
the majority of the ensemble. However, the oscillations of the pion 
fields are well within range, which is due to the rather small kinetic 
energy corresponding to (|20|). The bottom figures also show the 
ensemble evolution for the Gaussian initial condition (p^; mean 
values — crosses, 70% interval - dashed. The evolution is close to 
that from Vq = 3.3fm^. 


13 





















































































































